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Context. The International Gamma-Ray Astrophyiscs Laboratory (INTEGRAL) is discovering a large number of new hard X-ray 
sources, many of them being HMXBs. The identification and spectral characterization of their optical/infrared counterparts is a 
necessary step to undertake detailed study of these systems. In particular, the determination of the spectral type is crucial in the case 
of the new class of Supergiant Fast X-ray Transients (SFXTs), which show X-ray properties common to other objects. 
Aims. Our goal is to perform spectral analysis and classification of proposed counterparts to HMXBs in order to characterize the 
system they belong to. 

Methods. We used the ESO/NTT Sofl spectrograph to observe proposed IR counterparts to HMXBs, obtaining K s medium resolution 
spectra (R = 1320) with a S/N > 100. We classified them through comparison with published atlases. 

Results. We were able to spectrally classify the six sources. This allowed us to ascribe one of them to the new class of SFXTs and 
confirm the membership of two sources to this class. We confirmed the spectral classification, derived from optical spectroscopy, of a 
known system, 4U 1907-09, showing for the first time its infrared spectrum. The spectral classification was also used to estimate the 
distance of the sources. We compared the extinction as derived from X-ray data with effective interstellar extinction obtained from 
our data, discussing the absorption component due to the circumstellar environment, which we observed in four systems; in particular, 
intrinsic absorption seems to emerge as a typical feature of the entire class of SFXTs. 

Key words. X-rays: binaries - stars: supergiants - accretion - infrared: stars 
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High Mass X-Ray Binaries (HMXBs) are systems composed of 
an early-type massive star and an accreting compact object. All 
sub-groups of HMXBs involve OB type stars and are commonly 
found in the galactic plane and in the Magellanic Clouds, among 
their OB progenitors. 

The majority of the known systems are Be/X-ray Binaries 
(BeXRBs), consisting of a neutron star accreting matter from 
the circumstellar equatorial disc of a Be star. Most of them are 
transient, exhibiting short and bright outbursts (L x ~ 10 36 - 10 37 
erg s _1 in the case of Type I outbursts, generally close to the 
periastron passage of the neutron star; Lx > 10 37 erg s in the 
case of Type II outbursts). 

In the second major class of HMXBs, the Supergiant X-ray 
Binaries (SGXRBs), the counterpart is an early supergiant star, 
feeding the compact object with its radially outflowing stellar 
wind. As a result, the SGXRBs are, generally, persistent systems 
(L x ~ 10 36 erg s" 1 ). 

The five-year INTEGRAL data possibly reveal a different 
scenario. A recent subgroup has been proposed by Negueruela 
et al. (2006 1, named Supergiant Fast X-ray Transients (SFXTs): 
these objects, associated with a supergiant companion, occasion- 
ally undergo a short period of X-ray activity lasting less then a 
day, typically a few hours ( |Sguera et al.|2005) , with a very differ- 



ent behavior from those observed in other X-ray binaries. These 
outbursts show very sharp rises, reaching the peak of the flare 
in <1 hour. The decay is generally of a complex kind, with two 
or three further flares. The physical reason for these fast out- 
bursts is still unknown, although theoretical speculations would 
connect them to some sort of discrete mass ejection from the 
supergiant do nor (|G olenetsk ii et al.||20 03 ) or to wind variabil- 
ity ( |in't Z and 2005J), or to the possible presence of a second, 
equatorial, wind component ( |Sidoli et al.|2007| l . Due to high in- 
terstellar absorption and to the transient nature of these sources, 
SFXTs are difficult to detect, and in most cases, the sources had 
not been detected by previous missions. To date, six objects have 
been firmly characterized as SFXTs, but many other systems are 
likely candidates, and their number has grown rapidly since the 
launch of INTEGRAL ( [Winkler et al.|2003| , so that they could 
actually constitute a major class of X-ray binaries. 

Up to now, the INTEGRAL survey of the Galactic Plane 
and central regions has revealed the existence of more than 
200 sources ( |Bird et al.|[2007| |Bodaghee et~aT1[2007] l in the 
energy range 20-100 keV, with a position accuracy of 2' - 3', 
depending on count rate, position in the FOV and exposure. 
A large fraction of the newly discovered sources are found to 
be heavily obscured, displaying much larger column densities 
{Njj > 10 23 cm' 2 ) than would be expected along the line 
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* Based on observations collected at the European Southern 
Observatory, Chile (Programme IDs: 077.D-0568 and 079.D-0562) 



of sight (see Kuulkers 2005}. These sources were missed by 
previous high-energy missions, whose onboard instruments 
were sensitive to a softer energy range. Moreover, optical 
counterparts to these obscured sources are poorly observable 
due to the high interstellar extinction, with Ay in excess of up to 
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Table 1. NTT journal of observations. We report in the fourth column the net accumulated exposure time. Column five gives the 
obtained signal-to-noise ratio. The references list in the last column relates to the identification of the optical/infrared counterpart. 



Source 


K mag 


Start time (UT) 


Exp. time (s) 


S/N 


IR Counterpart 


Reference 


IGR J 16207-5 129 


9.1 


2006-07-14 23:19 


600 


100 


2MASS J16204627-5 130060 


Tomsick et al. (2006) 


IGR J 16465-4507 


9.8 


2007-05-26 05:03 


240 


100 


2MASS J16463526-4507045 


Zurita and Walter (2004) 


IGR J 16479-45 14 


9.8 


2007-05-26 05:11 


240 


100 


2MASS J16480656-4512068 


Kennea et al. (2005) 


AX J1841. 0-0536 


8.9 


2006-07-14 03:57 


600 


180 


2MASS J18410043-0535465 


Halpern et al. (2004) 


4U 1907+097 


8.8 


2006-07-15 07:37 


600 


130 


2MASS J 19093804+0949473 


Schwartz et al. (1980) 


IGR J 19 140+0951 


7.1 


2006-07-14 04:46 


360 


130 


2MASS J19140422+0952577 


in't Zand et al. (2004) 



~ 20 mag. 

In this context, the recent availability of infrared spec- 
troscopy has emerged as a strong tool to characterize these 
systems and, together with high-energy data, reveal the HMXB 
sub-class they belong to. This results in the identification of 
the mass transfer process of the system, with information about 
the intrinsic physics of the X-ray binary. The need for low 
energy data is particularly urgent in the case of SFXTs, which 
show X-ray properties common to other objects (such as RS 
CVs binaries and Low Mass X-ray binaries) and thus crucially 
require the spectral classification of their counterpart in order to 
be properly discerned. 

In this paper we present spectral analysis and classification 
of six HMXBs discovered (or re-descovered) by INTEGRAL. 
The selected IGR sources are the following: IGR J 16207-5 129, 
IGR J16465^507, IGR J16479-4514, AX J1841. 0-0536 and 
IGR J 19 140+0951. We also included the well known system 4U 
1907+09 since the spectral classification of its counterpart has 
been a matter of debate in the past, and no infrared spectra have 
been published up to now. The first three sources are located 
in the direction of the Norma-arm tangent region, the fourth in 
the Scutum-arm tangent region, the fifth and the last one in the 
Sagittarius arm tangent. In the next section we describe the ob- 
servations and data reduction; in Section [3] we report the ob- 
tained spectra, analyze their features and propose a classifica- 
tion; we calculate the interstellar hydrogen column density and 
estimate the distance to each source; in section|4]we discuss our 
results, before concluding. 

Preliminary results of our data analysis for AX J 184 1.0- 
0536 a nd IGR J 19 140+0951 were published in |Nespoli et~aT 
( |20Q7] >. 

2. Observations and data analysis 

We selected proposed counterparts, choosing sources observed 
by X-ray missions like XMM or Chandra, which can produce a 
very small error circle and facilitate the detection of the counter- 
part. 

Data were obtained in visiting mode during two observing 
runs, in July 2006 and May 2007 respectively, at the European 
Southern Observatory (ESO). The employed instrument was the 
SofI spectrograph (Moo rwood et al.| [l998), on the 3.5m New 
Technology Telescope (NTT) at La Silla, Chile. Table [T] reports 
the observation log, including, for each spectrum, the retrieved 
signal-to-noise ratio (S/N). 

We used the long slit spectroscopy mode, at medium 
resolution (R = 1320) with a K s grism and 1" width slit. The 
instrument large field objective provided a FOV of 4.92' x 4.92'. 
The sky had thin cirri on 2006 July 14th and 15th, while it was 



generally clear on 2007 May 26th. Seeing averaged between 
0.9" and 1.2", with the exception of the observation of IGR 
J 16749^45 14 which was performed with a seeing of 1.6". 

In order to ensure accurate removal of atmospheric features 
from the spectra, we followed a strategy similar to that outlined 
by |Clark & Steele| ( |2000[ ). At the telescope, we observed an 
AO - A3 III-V standard star immediately before or after each tar- 
get and a G2-3 V star once per hour in order to obtain very small 
differences in airmass (differences between 0.01 and 0.04 air- 
masses were accomplished). To compute the telluric features in 
the region of the H 1 21 661 A (Brackett-y line, or Bry), which is 
the only non-telluric feature in the A-star spectra, we employed 
the observed G-star spectra divided by the solar spectrum 1 prop- 
erly degraded in resolution. The dispersion solution obtained for 
the SofI spectra was also applied and the spectra of the A star, G 
star and the solar one were aligned in wavelength space. A tel- 
luric spectrum for each scientific target was obtained by patching 
into the A-star spectrum the ratio between the G star and the so- 
lar spectrum in the Bry region (we selected the range 21 590 - 
21 739 A). 

Typical integration times for standard stars were between 3 
and 7 minutes. 

Data reduction was performed using the IRAF 2 package, fol- 
lowing the standard procedure. We first corrected for the inter- 
quadrant row cross-talk, a feature that affects the SofI detec- 
tor; we then applied sky subtraction; we employed dome flat- 
fields and extracted one dimensional spectra. Wavelength cali- 
bration was accomplished using Xenon and Neon lamp spectra. 
Spurious features, such as cosmic rays or bad pixels, were re- 
moved by interpolation, when necessary. The reduced spectra 
were normalized by dividing them by a fitted polynomial con- 
tinuum. We corrected for telluric absorption, dividing each sci- 
entific spectrum by its corresponding telluric spectrum, obtained 
as described above. A scale and a shift factor were applied to 
the telluric spectrum, to best correct for the airmass difference 
and the possible wavelength shift; the optimum values for these 
parameters were obtained using an iterative procedure that min- 
imizes the residual noise. 

3. Results 

In this section we present the results of spectral classification and 
analysis for each target. The field of NIR spectral classification 

1 We used the NSO/Kitt Peak FTS solar spectrum, produced by 
NSF/NOAO 

2 IRAF is distributed by the National Optical Astronomy 
Observatories which is operated by the Association of Universities for 
Research in Astronomy, Inc. under contract with the National Science 
Foundation 
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is still very young and the level of required S/N and resolution to 
perform a quantitative profile analysis are very high (R ~ 12 000 
and S/N > 250), especially for young massive stars, as shown 
by Hanso net al.| ( |200 5b). The difficulty of the analysis depends 
on the few lines available and the relatively large uncertainty in 
their strength, due to their intrinsic weakness; moreover, some 
significant spectral regions, specifically through the 20 580 He I 
and the Bry features, pose systematic complications because of 
the strong telluric absorption. 

Under this premise, our analysis will be qualitative, based on 
the comparison with available NIR spectral atlases (Hanson et al 
[1996"" 



|2005a[ ). According to our estimation, this approach im- 
plies that the resulting spectral classification is provided with an 
uncertainty of no more than one luminosity subclass. The spec- 
tral type is precise up to one subtype. 

Some of the identified features exhibit a 10-20 A displace- 
ment with respect to the nominal values, consistent with the in- 
strumental resolution. Greater displacements, up to a maximum 
of 29 A, are found corresponding to lines characterized by a 
complex profile, such as Bry, or placed in a region of strong 
telluric absorption, such as both Hi 20581 A and Bry. 

In the next section we will classify the spectra, showing the 
features we were able to identify, with the corresponding equiva- 
lent width. Although it was pointed out ( jHanson et al.|1996] > that 
equivalent widths may show variations between stars, we report 
them for completeness. The values we calculated are affected 
by an average 15% error. For each object, a figure displays the 
spectrum we obtained, together with some comparative spectra 
from the atlas of |Hanson et al. (1996]) in a separate box. Table [2] 
summarizes the line identifications for all the observed targets. 

3.1. Spectral analysis and classification 
3.1.1. IGRJ1 6207-51 29 

The source was discovered by IBIS/ISGRI in the first Galactic 
Plane Survey performed by INTEGRAL ( |Bird et al.|[2004) l, 
which measured a flux of 3.8 + 0.3 mCrab in the 20-40 keV 
range. Instead, only the upper limit of < 4 mCrab was ob- 
tained in the 40-100 keV band. Subsequent Chandra obser- 
vations allowed the identification of the optical/infrared coun- 
terpart. This was associated with USNO-B1.0 0384-056 0857 
= 2MASS J 16204627-5 130060 by |Tomsick eTal] ([2006> . The 
power-law spectral fit provided Nh = 3.7*J x 10 cm" 2 and 
photon index equal to indicating an intrinsically hard 

source. The lower limit to the stellar temperature was estimated 
to > 18 000 K, revealing the presence of a very hot, massive star. 
From the fit of optical/IR spectral energy distributions, the dis- 
tance was estimated to be between 3-10 kpc (3-9 kpc in the case 
of a supergiant classification). Subsequent optical spectroscopy 
from Mas etti et al.| ( |2006[ ) refined the distance estimate to 4.6 
kpc. 

|Negueruela & Schurch| ( |2007| l, from optical observations, 
constrained the spectral type to earlier than B 1 . 

Figure[T]shows the K s spectrum we obtained, with identified 
spectral features marked. 

The spectrum shows no metal lines (no N m or C iv), strong 
He i 20581 A emission, Hei absorption at 21 126 A and mod- 
erately strong Bry absorption. The atomic transitions observed 
are the typical of OB star spectra. The He i 20581 A line is a 
prominent feature in supergiant stars, so that it is considered 
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Fig.l. K. s spectrum for 2MASS 16204627-5130060, the in- 
frared counterpart of IGR J16207-5129. The positions of identi- 
fied spectral features are marked by solid lines. 



an important tracer of stars with extended atmospheres. It 
becomes weak or even disappears in main sequence stars and it 
is observed in emission in B type supergiants, whereas it is in 
absorption in O type supergiants (Han son et al.|1996 >. The Hei 
21 126 A line is present in late O - early B spectra. 



For what was outlined above, and by comparing the relative 
strengths of the identified lines with those of Hanson etal. (1996 
|2005a| l, we estimate the spectral type of the IGR J16207-5129 
counterpart to be B 1 la. This allows us to classify the system as 
a SGXRB, as also inferred by |Negueruela & S church] 0Ml\ . 

3.1.2. IGR J1 6465-4507 

The source was discovered by INTEGRAL during its only ob- 



served X-ray flare on 2004 September 6th-7fh (Lutovinov et al. 
2004[ ). Using data of XMM/Newton, |Lutovinov et al.| ( |2005j 
found pulsations with period P s = 228 s in the X-ray flux and 
high photoabsorption, with Nh ~ 7 x 10 23 cirT 2 . The single star 
falling in the XMM/Newton error circle was identified by Zurita 
& Walter] ( |2004[ ) as the counterpart to the X-ray source, and as- 
sociated with 2MASS J 164635 26-4507045. 

From an optical spectrum, Negueruel a et al.| (12006) classi- 
fied the source as a Bl lb supergiant, subsequently refining their 
classification to B0.5 lb (Ne gueruela et al.(2 007). They also es- 
timated the distance of the source to 8 kpc. 

Our IR data suggest a different spectral classification. 

Figure [2] shows the K s spectrum we obtained, with identified 
spectral features marked. 

The spectrum shows Hei 20581 A in quite strong absorp- 
tion; a feature at 20730 A is possibly identifiable with a faint 
blend of the Civ transitions at 20690, 20780 and 20830 A, 
which appear in the atlas from |Hanson et al] ( |1996| l only in the 
very high S/N spectra. We observe Hei 21 126 and Bry, both in 
strong absorption. 

We agree with Negueruela et al. ( 2007} in the classification 
of the counterpart as an early supergiant. However, we suggest 
that the observed features point to an earlier type than Bl: in 
fact, the He 1 20 58 1 A line is seen in emission in B 1 supergiants, 
while it is in absorption in late-O supergiant. Comparing the rel- 
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Fig. 2. K s spectrum for 2MASS J 16463526-4507045, the in- Fig- 3. K s spectrum for 2MASS J 16480656-45 12068, the in- 
frared counterpart of IGR J16465-4507. The positions of identi- &ared counterpart of IGR J16479^514. The positions of iden- 



fied spectral features are marked by solid lines. 



tified spectral features are marked by solid lines. 



ative strengths of the outlined features with those from the at- 
lases of Hanson et al. ( 1996 2005a ), we refine the spectral clas- 
sification of the source to 09.5 la. We thus confirm the super- 
giant nature of the companion, which, together with the X-ray 
behavior of the system, classifies it as an SFXT, as proposed by 
Negueruela et al.|(|2006"ll. 



3.1.3. IGR J1 6479-451 4 



The source was discovered by INTEGRAL (Molkov et al. 2003 ) 
during an outburst. The X-ray spectrum is fitted with a power 
law with a high-energy cut-off, with spectral in dex T - 1 .4, and 
the column density is Nh = 12 x 10 22 cm -2 (Lutovinov et al 
[2005] l. 

IGR J 16479-45 14 showed short outbursts with very fast 
observed in 2003 by |Sguera et all (|2003]l. Its X-ray 



rises, 



behavior is thus typical of SFXTs (Negueruela et al. 2006), 



however the lack of an optical/infrared spectral classification 
prevented the possibility of enrolling it in this new class of 
objects. 
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Fig. 4. K. s spectrum for 2MASS 18410043-0535465, the infrared 
counterpart of AX J1841. 0-0536. The positions of identified 
spectral features are marked by solid lines. 



The counterpart to the source was identified by |Kennea et al.| 
( [2005] ) through SWIFT observations. 



Figure [3] shows the K s spectrum we obtained, with identified 
spectral features marked. 

The spectrum shows He 1 20 58 1 A in strong absorption, well 
recognizable although affected by a clear telluric residual. The 
difference we obtained in airmasses between the telluric standard 
star and the target was in this case very low, equal to 0.004: we 
thus suppose that the poor correction of telluric absorption is due 
to the passage of a cirrus during the observation of the telluric 
standard. 

We can also detect absorption at H ei 21 126, a weak Nm 
21 155 A emission (which, according to Hanson et al. (2005ai, 



could alternatively be C m), and moderately strong Bry absorp- 
tion. 

We conclude that the spectrum shows the typical features of a 
late O supergiant, especially the presence of He 1 21 126 (typical 
of late-O and early-B stars) in combination with He i 20 58 1 A, 
present in supergiant stars, and seen in absorption in late-0 su- 



pergiants. Through the comparison with the atlases from |Hanson 
et al. ( 1996 , 2005a ), we estimate the spectral type to be 09.5 lab. 
This result, together with those from X-ray data, allows us to af- 
firm that the object belongs to the new class of SFXTs. 

3.1.4. AX J1 841 .0-0536 

AX J1841. 0-0536 was discovered as a violently variable tran- 
sient by ASCA in April 1994 ( |Bamba et aLl[200T] ). The source 
showed multi-peaked flares with a sharp rise. Analysis of the 
ASCA data revealed that the source is a pulsar with P sp m - 4.7 s. 
The spectral fit provided a value of Nh - 3.2 x 10 22 cm -2 
( [Bamba et al.||2003[). A fast outburst observ ed by INTEGRAL 
was attributed by Halpern & Gottheff| ( |2004| ) to this source. 



A Chandra observation of the field revealed the counterpart 
to be 2MASS 18 410043-0535465 a re ddened star with weak 
Har in emission ([Halpern et al. 2004[), suggesting it was a 



Be star. Negueruela et al. (2006i, from optical spectroscopy, 



proposed the star is instead a luminous B0-1 type, although with 
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Fig.5. K. s spectrum for 2MASS J 19093804+0949473, the in- Fig. 6. K. s spectrum for 2MASS J19140422+0952577, the in- 
frared counterpart of 4U 1907+09. The positions of identified frared counterpart of IGR J19140+0591. The positions of iden- 



spectral features are marked by solid lines. 



tified spectral features are marked by solid lines. 



some uncertainty, classifying the system as an SFXT 

Figure |4] shows the K s spectrum we obtained, with identified 
spectral features marked. 

The spectrum shows He i 20581 A emission, accompanied 
by a spurious feature, possibly due to poor telluric component 
removal; we observe absorption at Hei 21 126, a weak Nm 
(Cm) 21 155 A emission; moreover, there is moderately strong 
Biy absorption. The side features of the Bry absorption profile 
are probably due to poor telluric correction, but they do not 
prevent us from measuring the equivalent width. 

The observed transitions are typical of an early supergiant, 



and by comparison with the atlases from Hanson et al. (1996 
2005a), we can conclude that the star is of Bl lb type. Together 



with X-ray properties, this NIR spectral classification allows us 
to confirm the nature of the system as an SFXT. 

3.1.5. 4U 1907+09 



The wind-accreting system 4U 1907+09 ( |Giacconi et al.|197T| is 
a known HMXB consisting of a neutron star in an eccentric (e 
0.28 ) 8.3753 day orbit around its companion, which has been op- 
tically identified as a highly reddened star ( |Schwartz et al.|1 980). 
The spectral classification of the counterpart to 4U 1907+09 has 
been matter of debate. The presence of X-ray flaring seen twice 
per neutron s tar orbit (|Marshall & Ricketts| 1 980|> had led some 
authors (e.g. |Makishima et al.||1984| |Cook & Page|[T987t [Tye 



1986| l to the hypothesis of a Be star companion. However, this 



classification would require a distance of <1.5 kpc, which is in 
contradiction to the significant interstellar extinction measured 
in optical observations by van Kerkwijk et al. ([1989 ), who also 
classified the counterpart as a B supergiant. Using interstellar 
atomic lines of Na I and K I, |Cox et aL] ( |2005] l set a lower limit 
of 5 kpc for the distance and proposed that the stellar companion 
is instead a 08-09 la supergiant with an effective temperature of 
30500 K, a radius of 26 Rq, a luminosity of 5 x 10 s Lq, and a 
mass loss rate of 7 x 10~ 6 M0 yr~' . 

Similarly to other accreting neutron stars, the X-ray contin- 
uum of 4U 1907+09 can be described by a power-law spectrum 
with an exponential turnover at 13 keV. The spectrum is mod- 



ified by strong photoelectric abs orption with a colum n density 
N H = 1.5 - 5.7 x 10 22 cm' 2 (e.g. |Cook & Page|l987) . 

We show for the first time an infrared spectrum of the source, 
which permits us to confirm the spectral classification as esti- 
mated from optical data. 

Figure [5] presents the K s spectrum we obtained, with identi- 
fied spectral features marked. 

The spectrum shows Hei absorption both at 20580 A and 
at 21 126 A, a weak N III (or C III) emission line at 21 155 A 
and strong Bry absorption (EW < 4 A), the typical features of 
an early supergiant. The presence of Hei 20580 in absorption 
strongly constrains the spectral type to a late O star. By com- 
parison with the atlases from Hanson et al. ( 1996 2005a 1, we 
conclude that the star is an 09.5 lab. We thus confirm and refine 
the previous spectral classification. 

3.1.6. IGR J19140+0951 

The INTEGRAL discovery of this source was reported by 
Hannikainen et al."| (2003 i. Observations with the Rossi X-Ray 
Timing Explorer (RXTE) revealed a rather hard spectrum, fit- 
ted with a power law of photon index 1 .6 and an absorption 
column density of Nh — 6 x 10 22 ctrT 2 ( |Swank & Markwardt 



2003). Timing analysis of the RXTE All- Sky Monitor (ASM ) 
data showed an X-ray period of 13.55 days ( Corbet et al.|2 004). 



Hannikainen et al. (2004) presented high energy spectral 



analysis of the period of the discovery, concluding that the 
source manifests two distinct spectral behaviors, the first 
showing a thermal component in the soft X-ray and hard X-ray 
tail, the second being harder and possibly originating from 
thermal Comptonization. This second, low-luminosity, behavior 
was confirmed to be the preferred state of the source (Rodriguez 
|et al|2005| ). 

The optical/inf rared counterpart to IG R J 19 140+0951 
was identified by in't Zand et al. ( 2006| ), from Chandra 
accurate position determination, as the heavily reddened 
2MASS 19140422+0952577. 

Figure [6] shows the K s spectrum we obtained, with identi- 
fied spectral features marked. The spectrum shows He 1 20 58 1 A 
emission, absorption at He 1 21 126 A, a weak N m (or C in) emis- 
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Table 2. A^-band line identifications and spectral classification: for each target spectrum, the observed wavelength, the corresponding EW and our 
proposed classification are reported. Positive EWs refer to features in absorption. 



Object 


Hei(2s 1 S-2p'P°) 
20581 A 


Civ(3p 2 P°-3d 2 D) 
20690, 20780, 20 830 A 


Hei(3p 3 P°- 
+ Hel (3p'P° 


4s 3 S), 21 120 A 
-4s'S), 21 130 A 


Nm (Cm) 
21 155 A 


Bry 
21661 A 




Spectral 
classification 


Wavelength 
[A] 


EW 
[A] 


Wavelength 
[A] 


EW 

[A] 


Wavelength 

[A] 


EW 

[A] 


Wavelength 

[A] 


EW 

[A] 


Wavelength 

[A] 


EW 

[A] 


ICR J 16207-5 129 


20595 


-2.90 






21 134 


1.23 






21642 


3.52 


Bl la 


IGR J 16465^1507 


20589 


2.72 


20730 


-0.42 


21 128 


2.31 






21 654 


3.25 


09.5 la* 


ICR J 16479^15 14 


20570 


3.10 






21 123 


1.23 


21 162 


-0.39 


21 640 


1.39 


09.5 Iab« 


AX J1841. 0-0536 


20599 


-0.94 






21 135 


1.64 


21 171 


-0.23 


21 652 


1.89 


Bl lb* 


4U 1907+09 


20602 


1.12 






21 130 


1.84 


21 170 


-0.25 


21641 


3.19 


09.5 lab 


IGRJ19140+0951 


20591 


-1.38 






21 136 


1.56 


21 167 


-0.33 


21632 


1.86 


Bl lab 



* Systems classified through this work as new or confirmed SFXTs. 



sion feature at 21 155 A, and moderately strong Bry absorption, 
typical features of an early supergiant. By comparison with the 
atlases from Hanso net al.| ( |1996| |2005a| l, we can conclude that 
the star is a Bl lab type. Together with X-ray properties, this al- 
lows us to confirm the nature of the system as an SGXRB. 



Our preliminary results were published in Nespoli et al. (2007), 
and recently confirmed by |Hannikainen et aL ^ 2007[ l, who, from 
K- and H-band spectra, constrained the spectral type to a B0. 5 
supergiant. They also estimated the distance of the source as 5 
kpc. 

3.2. Reddening and distance estimation 

From the identified spectral types, we obtained the intrinsic col- 
ors (J - K)o from Wegner (1994); we then calculated, from 
2MASS photometry, the instrumental color s (J-K)7mass, p rop- 
erly transformed through the formula from Carpenter ( 200 1 \ 3 to 
the Bessell & Brett ( 1988[ ) homogenized photometric system in 
order to estimate the infrared color excess E( J - K). Assuming 
the mean extinction law (R v = 3.1), from A V /E(J - K) = 
5.82±0.10 ( |Rieke & Lebofsky|1985] ), we obtained the total mea- 
sured visual extinction Ay and the corresponding hydrogen col- 
umn density value from Nh /Ay = 1 .79 + 0.03 x 10 21 atoms cm -2 
mag (Predehl & Schmitt 1995 ). We were thus able to compare 
the retrieved interstellar value of Nh with the one provided by 
X-ray data. In our calculation, we estimated errors through error 
propagation. Errors in the final values of Nh are mainly due to 
errors in the infrared colors and of the transformation between 
the two photometric systems. 

We also estimated the distance of the six sources, ap- 
plying the relation Mk = K + 5 - 5 log d - Ak- For each 
source, M v was obtained from our proposed spectral type 
( Wegner 2006 ), the intrinsic color index (K - V)q from Wegner 
(1994) was used in order to calculate M K , and the 2MASS K 
magnitude was employed. We derived A k from the r elation 
A A /E(J -K) = 2.4(A)- 1 - 75 , for A = 2.2 /mi ( |Draine|l989] ). 



The results of our calculations are given in Table [3] together 
with some crucial quantities used in the calculations or displayed 
for comparison. The distance estimation is mainly affected by 
the uncertainty in the value of the absolute magnitude My, which 
is due to two contributions, the errors given in the tabulated val- 
ues of My and the uncertainty in the spectral classification, from 
which the absolute magnitude is determined. The largest role is 
played by the errors in the mean tabulated values themselves (see 



3 In its updated version at |http://www.astro.caltech edu~jmc/2 mass/| 
/v3/transformafions/ 



Wegner 2006 for more details). The retrieved values for d must 
thus be assumed with prudence. 



4. Discussion 

Using near-infrared spectroscopy of six high-energy sources, 
IGR J16207-5129, IGR J16465^1507, IGR J16479-4514, 
AX J1841.0-0536, 4U 1907+097 and IGR J19140+0951, we 
classified their counterparts through comparison with published 
atlases. We found that all the observed systems have a supergiant 
companion. Our results, combined with information from X-ray 
data, is able for the first time to firmly include one source, 
IGR J16479-4514, in the newly discovered class of the SFXTs. 
Moreover, we can confirm with infrared data the identification 
of IGR J16465-4 507 and AX J1841.0-05 36 as SFXTs, as re- 
cently proposed by Negueruela et al. (2007 ) from optical spectra. 

From our spectral classification, we estimated the distance 
of the six sources. The retrieved values are consistent with the 
location of the sources in the Norma (IGR J16207-5129, IGR 
J16465-4507, IGR J16479^1514), Scutum (AX J1841.0-0536) 
and Sagittarius arm (4U 1907+097 and IGR J19140+0951) 
regions respectively. This determination, although affected by 
some uncertainty, can be considered an a posteriori confirma- 
tion of the proposed spectral classification. 

This work allowed us to calculate the extinction from IR 
data for the six systems. Recently, it has been pointed out (see 
Kuulkers||2005| |Chaty||2007] ) that INTEGRAL is revealing two 
new classes of supergiant HMXBs, the highly obscured HMXBs 
and the SFXTs. The first ones are characterized by strong intrin- 
sic absorption, the second by strong and short X-ray outbursts. 
High, variable, hydrogen column densities have in some cases 
been m easured for SFXTs (e.g. IGR J 11215- 5952: |Smith et al. 
2006al IGR J 1 739 1-302 1 : |Smith et al.|2006b| ), suggesting a pos- 
sible intrinsic absorption for this class as well, and marking a 
potential overlap between the two new classes. The origin and 
position (around the compact object only, or enveloping the en- 
tire system) of the absorbing material are still a matter of debate, 
and only multiwavelength studies are able to address the prob- 
lem, distinguishing between the absorption in X-ray and in the 
IR/optical bands. 

We calculated the effective interstellar extinction Ay and 
converted it into hydrogen column density Nh- Our results can 
be compared with the values obtained from X-ray data. If the two 
retrieved values are compatible within the corresponding errors, 
we face two possible scenarios: either the source of absorption 
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Table 3. For each observed source in the first column, the obtained spectral classification, intrinsic infrared colors, 2MASS pho- 
tometry, calculated infrared excess, hydrogen column density obtained from X-ray published measurements, effective interstellar 
column density obtained from our work and distance estimation are reported. See text for corresponding references. 



Source 


Spectral type 


(J ~ K) 


(J ~ KhMASS 


E(J - K) 


Nh from X-ray data 


Interstellar Nh 


Distance 






[mag] 


[mag] 


[mag] 


[10 22 cm" 2 ] 


[10 22 cm" 2 ] 


[kpc] 


IGR J 16207-5 129 


Bl la 


-0.12 


1.31 


1.43 


3 7 +L4 


1.53 ± 1.02 


6.1 (-3.5, +8.9) 


IGR J 16465^1507 


09.5 la 


-0.15 


0.69 


0.84 


72 ±6 


0.87 ±0.56 


9.5 (-5.7, +14.1) 


IGR J 16479-45 14 


09.5 lab 


-0.14 


3.19 


3.33 


12 ±4 


3.47 ±2.16 


2.8 (1.7, +4.9) 


AX J1841. 0-0536 


Bl lb 


-0.13 


0.80 


0.93 


3.2 


0.97 ± 0.64 


3.2 (-1.5, +2.0) 


4U 1907+097 


08.5 lab 


-0.17 


1.22 


1.39 


1.7-5.7 


1.45 ±0.87 


2.8 (-1.8, +5.0) 


IGR J 19 140+0951 


B0.5 lab 


-0.12 


1.50 


1.62 


~ 6° 


1.68+1.5 


1.1 (-0.8, +2.3) 



a A maximum value of 10.1 ± 0.2 x 10 22 cm 2 was reported by Rodriguez et al. 2005. 



is just the interstellar medium or there is a contribution from 
an extensive envelope around the whole binary system, if the 
extinction shows an excess of some orders of magnitude with 
respect to the estimated IS value. Conversely, if the reddening 
(measured from the IR colors) is low for the measured Nh from 
X-rays, this reveals the presence of an additional source of ex- 
tinction, which only affects the compact object in which the X- 
ray emission originates, and that can be regarded as the presence 
of absorbing material around it. 

In our case, the comparison reveals that for four systems, 
IGR J16465^507, IGR J16479-4514, AX J1841. 0-0536 and 
IGR J19140+0951, the extinction measured at high energy 
shows an excess of one or two orders of magnitude with respect 
to that obtained from IR data, probing that the material absorbing 
in the X-ray is concentrated around the neutron star. Our result 
is consistent with the detection of the Fe fluorescence line at 6.4 
KeV in the X-ray spectra of all the four sources ( ] Walter et al.| 
2006} |Bamba etal.|2001| |Rodriguez et al.|2005| l, which is con- 
sidered a signature of the dense spherical envelope around the 
compact object. 

Usually, the so called highly absorbed IGR sources are 
identified from their measured Nh being > 10 23 cirr 2 , i.e. one to 
two orders of m agnitude higher than the assumed Galactic value 
of ~ 10 22 cm 2 ( |Kuulkers|2005| l. The comparison of X-ray data, 
which are sensitive to the absorption from the environment of 
the compact object, with infrared data, which are in general 
only absorbed by the interstellar medium, can be a strong and 
alternative criterion to recognize this class of highly absorbed 
sources. 

In particular, all the three identified SFXTs show strong intrinsic 
absorption from circumstellar material, reported as well for 



other members of this new class (e.g. Negueruela et al.[2006 



Pellizza et al. 2006 1, and so emerging as a common feature 
of the group. This scenario is relevant to discussions about 
the physical mechanism, still unknown, which powers the 
fast outbursts that these sources undergo. Together with their 
X-ray fast transient nature, the high intrinsic absorption can 
be regarded as one of the reasons why this class has remained 
undiscovered until only recently. 

A more complete, multiwavelength investigation is neces- 
sary to reveal the nature of these systems, especially if it is able 
to find and physically describe correlations between their X-ray 
behavior and IR properties. In the case of SFXTs this is possibly 
the only way to explore the physical mechanism driving the fast 
outbursts. 



5. Conclusions 

From near-infrared spectroscopy of the six HMXBs we have 
found that: 

- the proposed optical counterparts were confirmed and the 
spectral classification of the sources provided; 

- one source, IGR J16479^+514, was added to the SFXTs 
and the confirmation of IGR J 16465-4507 and AX J1841.0- 
0536 as members of the class was proven with infrared data; 

- the comparison between Nh obtained from X-ray data and 
interstellar extinction from our data showed for four systems 
(IGR J16465-4507, IGR J16479-4514, AX J1841.0-0536 
and IGR J19140+0951) the presence of an absorbing enve- 
lope, strictly confined around the compact object; 

- all the three identified SFXTs are intrinsically absorbed, sug- 
gesting that this might be a characteristic of the class; 

- the distance estimation, compatible with the location of the 
sources in the respective galactic arms, is a possible confir- 
mation of the spectral classification provided here. 
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